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Abstract For the two sixth-order polynomials a(s) and 6(s), Hurwitz stability of their convex 
combination is necessary and sufficient for the existence of a polynomial c(s) such that c(s)/a(s) and 
c{s)/b{s) are both strictly positive real. Our reasoning method is constructive, and is insightful and 
helpful in solving the general robust strictly positive real synthesis problem. 
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The strict positive realness (SPR) of transfer functions is an important performance specification, 
and plays a critical role in various fields such as absolute stability/hyperstability theory [12, 17], 
passivity analysis [9], quadratic optimal control [2] and adaptive system theory [13]. In recent years, 
stimulated by the parametrization method in robust stability analysis [3, 5], the study of robust strictly 
positive real systems has received much attention, and great progress has beem made [1, 4, 6, 7, 8, 
10, 11, 14, 15, 16, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30]. However, most results belong to 
the category of robust SPR analysis. Valuable results in robust SPR synthesis are rare. The following 
fundamental problem is still open [1, 4, 8, 10, 11, 14, 15, 16, 18, 22, 24, 26, 27, 28, 29]: 

Suppose a{s) and b(s) are two n-th order Hurwitz polynomials, does there exist, and how to find a 
(fixed) polynomial c(s) such that c(s)/a{s) and c{s)/b{s) are both SPR? 

By the definition of SPR, it is easy to know that the Hurwitz stability of the convex combination 
of a(s) and b{s) is necessary for the existence of polynomial c(s) such that c(s)/a(s) and c{s)/b{s) 
are both SPR. In [10, 11, 16], it was proved that, if a{s) and b{s) have the same even (or odd) parts, 
such a polynomial c(s) always exists; In ]1, 10, 11, 14, 19, 20, 22, 23, 26, 27], it was proved that, 
if n < 4 and a(s),6(s) G K {K is a stable interval polynomial set), such a polynomial c(s) always 
exists; Recent results show that [19, 20, 22, 23, 24, 26, 28, 29], if n < 5 and a(s) and 6(s) are the 
two endpoints of the convex combination of stable polynomials, such a polynomial c(s) always exists. 
Some sufficient condition for robust SPR synthesis are presented in ]1, 4, 8, 14, 19, 20, 23], especially, 
the design method in [19, 20] is numerically efficient for high-order polynomial segments and interval 
polynomials, and the derived conditions are necessary and sufficient for low-order polynomial segments 
and interval polynomials. 

This paper shows that, for the two sixth-order polynomials a(s) and b(s), Hurwitz stability of 
their convex combination is necessary and sufficient for the existence of a polynomial c(s) such that 
c(s)/a(s) and c{s)/b{s) are both SPR. This also shows that the conditions given in [19, 20] are also 
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necessary and sufficient, and the open problem above has a positive answer for the case of sixth-order 
polynomial segment. Our reasoning method is constructive, and is useful in solving the general robust 
SPR synthesis problem. 

In this paper, P" stands for the set of n-th order polynomials of s with real coefficients, R stands 
for the field of real numbers, d{p) stands for the order of polynomial p(-), and P" C P" stands for 

the set of n-th order Hurwitz stable polynomials with real coefficients. 

In the following definition, p(-) G P™, q{-) G P", /(s) = p{s)/q{s) is a rational function. 

Definition 1 [25] /(s) is said to be strictly positive real(SPR), if 

(i) d{p)=diqy, 

(ii) /(s) is analytic in Re[s] > 0, (namely, q{-) G i?" ); 

(hi) Re[/(ju;)] >0, Vw G R. 

If f{s) = p{s)/q{s) is proper, it is easy to get the following property: 

Property 1 [7] If /(s) = p{s)/q{s) is a proper rational function, q{s) G i7", and Vw G 
P, Re[f {ju)] > 0, then p{s) G i7” U 

The following theorem is the main result of this paper: 

Theorem 1 Suppose a(s) = s®-|-ais^-|-a 2 s‘*-l-a 3 s^-|-a 4 s^-|-a 5 s-|-a 6 GH^,b{s) = + 

bss^ b^s^ -\- b^s -I- feg S the necessary and sufficient condition for the existence of a polynomial 

c(s) such that c(s)/a(s) and c(s)/6(s) are both Strictly Positive Real is 

Xb{s) + {1 - X)a{s) g77®,Ag [0,1]. 

Since SPR transfer functions enjoy convexity property, by Property I, we can easily get the nec¬ 
essary part of the theorem. 

To prove sufficiency, we must first introduce some lemmas. 

Lemma 1 Suppose s® -I- -|- 025 "* -I- 035 ^ -I- 045 ^ -I- a^s -|- Og G then the following quadratic 

curve is an ellipse in the first quadrant of the x-y-z-p space: 

{ {a 2 X z — aiy — 03 )^ — 4(ai — a;)(a 5 -I- asy aip — a 2 Z — 04 ^) = 0 
qqx a4Z — asp — a^y = 0 
osP - CLqZ = 0 

and this ellipse is tangent with the line 


qqx a^z — a^p — a^y = 0 
a^p — uqz = 0 
oi — X = 0 


at 


X = ai, 

QSQ-qQ-i — Q-5 0 , 40,2 Q-1~I~Q5Q4Q-3~I~Q3 0-6^2 0-1 — O-^O-Q 


y = 

Z = -05 


a^a20,i—aQa^—a^a3 
— 05(1401-t-flg-t-a3a6 fli 


P = 


O 5 O 2 O 1 — OgO,, — O 5 O 3 
® — 05040i-t-0g-t-030g0i ' 


tangent with the line 


qqx 042 — 03P — a^y = 0 
a^p — gqz = 0 

05 + 032/ + aiP ~ o,2Z — a^x = 0 
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at 


— Cl^-\-Cl^Cl^(ll “l“Cl 3 ClQ~\~Cl^ 0-2 ®3 — < 3,5 Gg <3-4 —2t35t33(36(3]^ 


a?ao —<3ct34+a5a? ai — 2a5a6aia2 — <33<36ait34+a?a?+a5a3a6 — a5a3a4a2+<3q<36t32 ’ 
^ ^ ^ 22 2z^2 ^ 

_ a6a4aq — a 3 aQai — 03050 ^ —a^aQ-\-a^a4a2 

- 9^*^ “5 1 

a|a|—a5a4+a5a|ai—2a5a6aia2 —<33a6(3i(34+t3ga^+a5a3a6 —<35a3a4a2+<33t36a2 ’ 


<32<3g — CLqCLqCLx — <35 <34 <33 +< 330-6 


^ 0^01 — 0^04 +050^01 — 205050102 —0306Oi 04 + 0^0^+050305—05030402+o|0602 ’ 

02 0 ? — O5O6O1 —O5 O4O3+O? 06 

JD (Xf\ ^ ^ ^ • 

^ 0^02—0^04+050^01 — 2O5O6O1O2 —03060104+0^0^+050306 —O5O3 04 02+0^0602 


Proof Since a(s) is Hurwitz stable, Lemma 1 is proved by a direct calculation. 

Lemma 2 Suppose s® + Ois^ + 025 "* + 035 ^ + 045 ^ + a^s + 036 then the following quadratic 
curve is an ellipse in the first quadrant of the x-y-z-p space: 


(05 + a^y + aip — a 2 Z — 04 ^)^ — 4 (a 2 X + z — aiy — a^)(aQX + 04 ^ — a^p — a^y) = 0 

oi — X = 0 

a^p — qqZ = 0 


and this ellipse is tangent with the line 


oi — X = 0 

a^p — uqz = 0 

a 2 X + z — aiy — 03 = 0 


at 


X = ai, 

_ — O5+O5O4O1 —O5O2O1 +O5 0203+050 ]^02 —Oi 06 O3 

y —020105+0^06+0305 - 

_ O3 —03 02 01 +04 of —0105 


z = 


-020105+0^05+0305 


p = 


Og —O3O2O1 +Q 4 0 -|^ — Ol O5 


4(120105+0^06+0305 


tangent with the line 


at 


Ol — X = 0 

osp — oez = 0 

oex + 042 — osp — a^y = 0 


X = tti, 

0^0501 — 03O4O105 — 0^04+030505 — O2O1 O5 05+O4 Og 

—0^02+050403—0^06+050601 ' 


y = 


05O4O1 —010503 —Og _ 

—0^02+050403—0^06+050601 ’ 


p = ae 


O5O4O1 —O1O5O3 — O5 
—0^02+050403—0^06+050601 ■■ 


Proof Since a(s) is Hurwitz stable, Lemma 2 is proved by a direct calculation. 

Lemma 3 Suppose s® + ois^ + 025 "* + 035 ^ + 045 ^ + a^s + 05 G then the following quadratic 
curve is an ellipse in the first quadrant of the x-y-z-p space: 

{ (oex + 04 ^ - a^p - osp)^ - 4(05 + a^y + oip — a2Z - aix){a^p — oez) = 0 

Ol — X = 0 

a 2 X + z - oip - 03 = 0 
and this ellipse is tangent with the line 
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at 


oi — a: = 0 

a2X + z — aiy — 03 = 0 

05 + 032/ + oip — a^z — a4X = 0 


X = ai, 

CL^CL^-\-CL^CL2^1'\'0-\0.Q — 0-2 <13 — CL^CL20-\ 
y ai—a:^a2CL\-\-CL4CLi—CLiCL^ ’ 

2 3 3 2 

(l2<llQ':^ — Cl2CliCl5~\~2o,i(l3Cl5-\-<22 dg — dg — dg d 4 d-|^ 

dg—d 3 d 2 di+d 4 d^ — di dg ’ 

2d5 d4di -|-d4 di d2 dg — dg — d5d2di-|-d5d2 dg — didgd3-|-dgd2d2 — d4dg — d^d^ 
^ dg—dg d2 d -^d4 d —dj^dg 


tangent with the line 


ai — X = 0 

a2X z — aiy — as = 0 
a^p — a^z = 0 


at 


X = ai, 

dgdgdi —d5d4d2d]^-l-d5 d4 d3-|-d3dgd2d]^ —d^dg 
^ dg d4di — didgdg —d^ ’ 

\ _ _ ^ —d2didg+d?dg+d3dg 

z — a s '2 ^ 

^ dgd4di —didgd3 —dg ' 

_ — d2didg+d^dg+d3dg 

P — a^ •} • 

V ^ dgd4di —didgdg —dg 

Proof Since a(s) is Hurwitz stable, Lemma 3 is proved by a direct calculation. 

For notational simplicity, denote 

figi := {{x,y,z,p)\ {a2X + z - aiy - 03)^ - 4 (ai - x){a5 + a^y + oip - a2Z 
—a^x) < 0, aex + a^z — asp — a^y = 0, a^p — aez = 0} 

^e 2 := {ix,y,z,p)\ (05 + asy + aip - a 2 Z - a^xY - 4 (a 2 X + z - aiy - as){aQX 

+a 4 Z — asP — a^y) < 0, oi — x = 0, a^p — agz = 0} 

^e 3 •= {( 2 ;,y, Z,p)\ {aex + a^z - a^p - a^yY - 4(05 + asy + aip - a 2 Z 

—a4x){a5p — uqz) < 0 , ai — X = 0 , a2X + z — aiy — 03 = 0 } 

^ei ■= {(x,y,z,p)l (b 2 X + z-biy- bsY - 4(6i - x )(&5 + h^y + hip - h 2 Z 
-h4x) < 0 , b^x + 64 Z - hsp - hsy = 0 , b^p - b^z = 0 } 

^e2 := {{x,yjZ,p)\ (65 + bsy + bip - b2Z - b 4 xY - ‘^{b2X + z - biy - bs){bQX 

+h 4 Z - bsp - bsy) < 0, 5 i - x = 0, b^p - b^z = 0} 

^e 3 := {{x,y,z,p)\ {bex + b4Z - bsp - b^yY - 4(&5 + b^y + bip - b 2 Z 

-&4x)(&5p — bez) < 0 , 61 - X = 0 , 62X + z - biy - 63 = 0 } 

In what follows, {A, B) stands for the set of points in the line segment connecting the point A and 
the point B in the x-y-z-p space, not including the endpoints A and B. Denote 

:= {(x,y,z,p)|(x,y,z,p) G {A,B) U {A,C) U {B,C),\/A G ^^i^B G L!:2,VC G 11:3} 

D" := {{x,y,z,p)\{x,y,z,p) G {A,B) U {A,C) U {B,C),\/A G G G ^^3} 

Lemma 4 Suppose a(s) = s® + ois® + 025"* + 035^ + 045^ + a^s + 05 G iL®, b{s) = s® + 6is® + 
b2s'^ + bss^ + &4S^ + 65S + &6 G if n Y 4 >^ take {x,y,z,p) G fl“ D and let c(s) := 
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c('S) c(s) 


(x — 5 ) 5 ^ H- ys^ + zs^ -\-ps-\- s is a sufficiently small positive number), then for —^ and , . . , 

a(5) b{s) 

we have Vw G R, 1 > 0 and Re[ 77 ^^] > 0. 

o(jw) b{juj) 

Proof Suppose {x,y,z,p) G Pi let c(s) := s® + (cc — e)s^ + ys^ + zs^ +ps + e,£ > 0,£ 
sufficiently small. 

Vw G R, consider 


Re[ 


c(jw) 1 


1 


[(oi — x)uj^^’ + {a 2 X + z — aiy - 03 ) 0 ;^ + (05 + asy + aip 


aUcoY |a(jo;)P' 

—a 2 Z — a4x)a;® + (ogx + 042: — asp — a 5 y)uj‘^ + {a^p — aQz)uP‘ 
-\-£{lo^^ — a 20 J^ + (04 — l)ui^ + (— clq + a 2 ')x>‘^ — a4W^ + cie)] 


In order to prove that Vw G i?, Re[d^—> 0, let t = uj'^, we only need to prove that, for any 

a{juj) 

£ > 0 ,£ sufficiently small, the following polynomial /i(t) satisfies 


fi{t) := t[(ai - xY'^ + (a2X + z - aiy - + (as + 03?/ + aip 

—a2Z — aix)t^ + (agx + 04^ — a^p — a^y)t + {a^p — aez)] 

-t-£(t^ — (12^^ (a4 — l)t^ (—clq 02)^^ — ci4.t ag)] > 0 , Vt G [ 0 , -t-00). 


Since (x, y, z,p) G fl“, by the definition of fl°', it is easy to know that 

gi{t) := (oi — x)t‘^ + {a2X + z — aiy — + (as + a^y + aip — a2Z — aix)t^ 

+ (a6X + 042; — a^p — a^y)t + [a^p — a^z) > 0 , Vt G [ 0 , +00). 


Moreover, we obviously have /i(0) > 0, and for any £ > 0, when t is a sufficiently large or sufficiently 
small positive number, we have fi(t) > 0 , namely, there exist 0 < ti < <2 such that, for all £ > 0 , 
t G [0,ti] U \t 2 , + 00 ), we have /i(t) > 0. 

Denote 

Ml = inf tgi{t), 
te[ti,t2] 

Ni = sup \Y — 02 ^^ + (04 — l)t^ + (—06 + a 2 )Y — ait + oe]. 


Then Mi > 0 and IVi > 0. Choosing 0 < £ < 


Ml 


by a direct calculation, we have 


fi{t) := t[{ai - x)t^ + {a 2 X + z - aiy - asY^ + (05 + asy + aip 
—a2Z — aixY'^ + {a^x + a^z — a^p — a^yY + (osP — a^z)] 

+£{t^ — a 2 t‘^ + (04 - l)t^ + (-06 + 02 )^^ - ait + oe) >0, Vt G [0, + 00 ). 


Namely, 

Vw G i?,Re[4^] > 0. 
a{juj) 

Similarly, since (x, y, z,p) G there exist 0 < fo < ^4 such that, for all £ > 0, t G [0, t 3 ]U[t 4 , + 00 ), 
we have / 2 (t) > 0,where 

/ 2 (t) := t[{bi - xY'^ + {h 2 X + z-hiy - b^t^ + (6s + hv + hp 
-b2Z - bixY^ + {bex + biz - b^p - b^yY + (6sP - b^z)] 

+£(t^ - 626"^ + (64 - l)t^ + {-b& + 62)6^ - bit + be)] 

Denote 

52(6) := (61 - x)t'^ + {b2X + z-biy- b^t^ + (6s + b^y + bip - 622; - bix)t^ 

+ {bex + 64 Z - bep - bey)t + {b^p - bez), 
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M 2 = inf tg 2 {t), 
ie[t3,t4] 

N 2 = sup \t^ - b 2 t'^ + (64 - + (-^6 + b 2 )t^ - bit + hn], 

te[t3,t4] 

M 2 

Then M 2 > 0 and N 2 > 0. Choosing 0 < e < , we have 

N 2 

\/uj € R,Re[^^] > 0 . 
b{juj) 

r ^2 T 

Thus, by choosing 0 < £ < minj-—, - 5 ^}, Lemma 4 is proved. 

Ni N 2 

Lemma 5 Suppose a(s) = s® + ois® + 025"^ + 035^ + 043^ + ass + ae G H^, b{s) = s® + 61 s® + 
&2s‘* + 635^ + bis'^ + b^s + be G if A6(s) + (1 — A)a(s) S i?®, A G [ 0 , 1 ], then 0 “ D 7^ (/) 

Proof If VA G [0, l],aA(s) := Xb{s) + (1 — A)a(s) G H^, by Lemmas 1 - 3 , for any A G [0,1], 
^ei; ^e 2 ^e 3 9'^’® three ellipses in the first quadrant of the x-y-z-p space, denote 

:= {{x,y,z,p)\{x,y,z,p) G {A,B) U iA,C) U (i?,C'),VAl G G G 


Apparently, when A changes continuously from 0 to 1, will change continuously from to 
■''^hl change continuously from 0 “^ to change continuously from ^‘^2 to fle 2 j s-od 

flg 3 will change continuously from Algg to flgg. 

Now assume 17“ n = </>, by the definitions of 17“ and 17^, and Lemmas 1-3, du > 0, u > 0, u 7 ^ 
oi, M 7 ^ &i, and 3k G {1, 2, 3}, such that the following plane in the x-y-z-p space 


I : 


X y z 

- 1 - 1 - 

U V w 



r 


separates 17“ and 17**, meanwhile, I is tangent with Og^, 17“2,17g3 and 17gj, simultaneously (or tangent 
with 17 ^ 3 , 17 ^ 2:^^63 and 17“^ simultaneously ). 

Without loss of generality, suppose that I is tangent with Og^, 17“2,17“3 and 17g^ simultaneously. 

Since I is tangent with 17“3,17“2 and 17“3 simultaneously, a(s) is Hurwitz stable and u > 0,u 7 ^ 

ai,u > 0 , by a lengthy calculation, we get that the necessary and sufficient condition for I being 

tangent with 17“^, 17“2 and 17“3 simultaneously is 

uv — aiv — a 2 uv -\- a^v -\- a^uv — a^v — agu = 0,w = —uv, r = —v (1) 

Since w = —uv, r = — I is tangent with 17g^, by a direct calculation, we have 

uv^ — biv^ — b2uv^ -\- b^v^ -\- biuv — b^v — b^u = 0 ( 2 ) 


^^From (1) and (2), we obviously have VA G [0,1], 

uv — a\iv — a\ 2 uv -\- ax^v -\- ax^uv — ax 5 V — axeu = 0,w = —uv, r = —v (3) 

where axi ■= ai-\-X{bi — ai),i = 1, 2,3,4, 5, 6 . (3) shows that I is also tangent with 17“^ (VA G [0,1]). But 
I separates 17“^ and 17g^, and when A changes continuously from 0 to 1, 17“^ will change continuously 
from 17“^ to 17g^, which is obviously impossible. This completes the proof. 

From Theorem 2.4 in [20], or the proof of Lemma 5 in [24], we have 

Lemma 6 Suppose a{s) = s® -|- ois® -I- 025 "^ -I- -\- 043 ^ -I- 055 -I- 05 G , b{s) = s® -I- 61 s® -I- 

& 2 s‘* -I- bzs^ biS^ b^s -Gbe G H^, c(s) = -I- xs'^ -\- ys^ -\- zs^ -\-ps-\-q, if Vw G R, Re[ *"^/.^| ] > 0 and 

a(jw) 

Re[ w'^^! ] > 0 , take 
b{juj) 
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c (s) := c(s) + 6 ■ d{s), 6 > 0, 6 sufficiently small 


c fsl c (5) 

(where d(s) is an arbitrarily given monic sixth-order polynomial), then —and -ttt' ^^e both 

a(s) b[s) 

strictly positive real. 

The sufficiency of Theorem 1 is now proved by combining Lemmas 1-6. 

Remark 1 From the proof of Theorem 1, we can see that this paper not only proves the existence, 
but also provides a design method. 

Remark 2 The method provided in this paper is constructive, and is insightful and helpful in 
solving the general robust SPR synthesis problem. This subject is currently under investigation. 

Remark 3 Our results can easily be generalized to discrete-time case. 

c(s) 


c( s ) c( s ) 

Remark 4 If and are both SPR, it is easy to know that VA G [0,1], 


is also SPR. 


a(s) b{s) 


Xa{s) + (1 — A)6(s) 


Remark 5 The stability of polynomial segment can be checked by many efficient methods, e.g., 
eigenvalue method, root locus method, value set method, etc. [3, 5]. 
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